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Abstract 

The brain of mammals are divided into different cortical areas that are 
anatomically connected forming larger networks which perform cognitive 
tasks. The cat cerebral cortex is composed of 65 areas organised into the vi¬ 
sual, auditory, somatosensory-motor and frontolimbic cognitive regions. We 
have built a network of networks, in which networks are connected among 
themselves according to the connections observed in the cat cortical areas 
aiming to study how inputs drive the synchronous behaviour in this cat 
brain-like network. We show that without external perturbations it is possi¬ 
ble to observe high level of bursting synchronisation between neurons within 
almost all areas, except for the auditory area. Bursting synchronisation ap¬ 
pears between neurons in the auditory region when an external perturbation 
is applied in another cognitive area. This is a clear evidence that pattern 
formation and collective behaviour in the brain might be a process mediated 
by other brain areas under stimulation. 
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1. Introduction 


The nervous system of mammals is responsible for collecting and process¬ 
ing information, where the signals are sent by neurons jili] . The propagation 
of neural signals occurs through electrical and chemical synapses, as a result 
of the difference in electric potential between the exterior and the interior of a 
neuron [2J. Neurons connect to each other forming complex layered structures 
[3j. The different cortical layers have particular distributions of neuronal cell 
types, as well as connections with other cortical and sub-cortical regions j4|. 
The mammalian brain is composed of distinct areas, the cerebellar cortex, 
and non-cortical nuclei. The cortex presents fundamental divisions such as 
the hippocampus formation, the olfactory cortex, and associated areas 0. 

In this work we consider the cat cerebral cortex. Scanned and collabo¬ 
rators 80 have relevant results related to the cortical system of the cat. 
They showed the connection organisation, and reported that there are 1139 
corticocortical connections among 65 cortical areas. The cortical areas are 
organised into four connectional clusters, corresponding to visual, auditory, 
somatosensory-motor, and frontolimbic areas js]. 

Here we focus on dynamical features such as bursting synchronisation 
and desynchronisation [9|. Bursting synchronisation are thought to play 
relevant roles in information binding in the mammalian brain [10]. However, 
bursting synchronisation may be associated with pathologies like seizures (ll| 
or Parkinson’s disease [l2[. For this reason, studies about synchronisation 
are of great interest to neuroscience. 

Our purpose in this work is to study the formation of patterns of synchro¬ 
nisation and desynchronisation in a neural network model of the cat brain, 
using the matrix of corticocortical connections in the cat 0. The matrix 
represents the densities of connections in 65 cortical areas is on undirected 
wei ghte d adjacency. We describe each cortical area as a small-world network 
13, 14, 13. Small-world networks have been proposed to be an efficient solu¬ 
tion for achieving phase synchronisation of bursting neurons 16]. In addition 


they have been found to be linked to different levels of models of the brain. 
Small-world networks have been intensively investigated in computational 


neuroscience 


17 . The characterisation can be made on two basic levels: a 


microscopic, neuroanatomic level, and a macroscopic, functional level. Stud¬ 
ies at the former level are limited to those few examples in which there is 
available data on the neuronal connectivity, as the worm C. Elegans, which 
is considered one of the simplest and most primitive organisms that shares 
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essential biological characteristics of the more complex species 18]. More 


over, there have been studies of large-scale anatomical connection patterns 
of the human cortex using cortical thickness measurements from magnetic 
resonance imaging [10. The human brain anatomical network at this level 
has an average path length and a clustering coefficient with values presented 
by networks with small-world property. Stam and collaborators, in a study 
of functional brain networks, observed that there is a loss of small-world net¬ 
work characteristics in patients with Alzheimer’s disease, in particular with 
an increase of the average path length with no significant changes in the 


clustering coefficient 20 


At the macroscopic level of description of neural networks, the use of non- 
invasive techniques as electroencephalography, functional magnetic resonance 
imaging and magnetoencephalography provides anatomical and functional 
connectivity patterns between different brain areas (21, 22]. This information 
provides a way to study the brain cortex, considering the latter as being 
divided into anatomic and functional areas, linked by axonal fibers. Scannell 
and coworkers have investigated the anatomical connectivity matrix of the 
visual cortex for the macaque monkey and the cat 0,0. In both cases the 
values of the average path length and clustering coefficient are in accordance 
with expected small-world properties 


Each node in our brain model is described by the Rulkov model 23 


a discrete time system with two dimensions. The low-dimensionality of this 
map allow us to study large neural networks of approximately 10,000 neurons. 
This model has been extensively tested and it reproduces well the main 
features of realistic neural models. It has been used in studies about control 
of bursting synchronisation [20, phase synchronisation in clustered networks 
[l0 . and suppression of bursting synchronisation [20 . 

Our main goal is to show that pattern formation and collective behaviour 
in the brain might be a process mediated by other brain areas under stimula¬ 
tion. We also show that external perturbations induce synchrony behaviour 
in cognitive areas of the cat cerebral cortex. Bursting synchronisation ap¬ 
pears between neurons in the auditory region when an external perturbation 
is applied in another cognitive area. We consider perturbations that activate 
neurons in accordance with experimental results in that pulses of blue light 
are capable to induce neuronal spikes (26[, 27 . 


This paper is organised as follows: in Section II we introduce the network 
of Rulkov neurons and the cat brain matrix. In Section III, we study the 
phase synchronisation of the cognitive areas according to electrical and chem- 
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ical synapses. In Section IV, we analyse the effect of an external perturbation 
on the synchronisation. In the last Section, we draw the conclusions. 


2. Network of Rulkov neurons 


There is a wide range of mathematical models used to describe neuronal 
In this work we consider the phenomenological model 


activity [23|, [281,1291, 30 
proposed by Rulkov 


a 

%n+l Z ' 2 ^ (f) 

1 + 4 

Vn+1 = Vn~ Cr(x n ~ p), ( 2 ) 

where x n and y n are the fast and slow dynamical variables, respectively. The 
parameter a affects the spiking time-scale, and we choose values in that the 
time series of x n presents an irregular sequence of spikes. The parameters a 
and p describe the slow time-scale. Figured] shows the time evolution of the 
fast and slow variables, where ilk is used to denote when neuronal bursting 
starts. 



Figure 1: Time evolution of the (a) fast and (b) slow variables in the Rulkov map. 
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Figure 2: (Colour online) Density of connections between cortical areas classified as absent 
of connection (white), sparse (red), intermediate (blue), and dense (green). 


We build a neuronal network considering the organisation of corticocor- 
tical connections obtained by Scannell and collaborators ( 3 , 0 ] - Figure [ 2 ] 

exhibits the matrix description of the corticocortical connectivity of the brain 
in accordance with Ref. [8[. In this work, the matrix elements have connec¬ 
tions weighted 0, 1, 2, or 3, where 0 represents the absence of a connection 
between two brain areas, 1 represents a sparse connection, 2 represents an 
intermediate connection, and 3 is dense connections. In Figure [2] the colours 
are representing the weights white for no connections, sparse connections in 
red, intermediate connections in blue, and dense connections in green. All 
areas are grouped into 4 cognitive regions: visual, auditory, somatosensory- 
motor, and frontolimbic. The visual region has 18 cortical areas, the auditory 
has 10 areas, the somatosensory-motor has 18 areas, and the frontolimbic has 
19 cortical areas. 

The cortical areas contain a characteristic distribution of neuronal cells 
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and connections. Supporting evidences that there are small world properties 
at different level models of the brain 3ll.l32| we consider each cortical area as a 
small world network. A small-world has an average distance among neurons 
like a random network, while the degree of clustering is comparable to a 
regular network [l4|. A small-world network has typically an average distance 
between sites comparable to the value it would take on for a random network, 
while retaining an appreciable degree of clustering, as in regular networks. 
Watts and Strogatz obtained small-world networks from an otherwise regular 
lattice with local connections, to which non local connections were added 
by randomly rewiring a small fraction of the local connections jlij]. An 
alternative procedure was proposed by Newman and Watts, who inserted 
randomly chosen shortcuts in a regular lattice, instead of re-wiring local 
links into non-local ones [33j. We build small-world networks according to 
the procedure proposed by Newman and Watts. Each small-world network 
has 100 neurons and %5 of shortcuts. The connections among the small-world 
networks obey the corticocortical connectivity of the cat (Fig. [2]) so that 2 
areas connected with weight equal to 1 (red) have 50 randomly connections, 
areas connected with weight equal to 2 (blue) have 100 randomly connections, 
and areas with weight 3 (green) have 150 connections. 

The coupling between neurons happens by means of electrical or chemi¬ 
cal synapses. The local connections between neurons within each small-world 
network are described by electrical synapses. The shortcut non-local connec¬ 
tions between neurons within each small-world network and the connections 
among cortical areas are described by chemical synapses. Chemical synapses 
may be excitatory or inhibitory. We consider that 75% are excitatory and 
25% are inhibitory f34| . 

The dynamic behaviour of the neuronal network with electrical and chem¬ 
ical connections is governed by the following equations 


x. 


ihv) 

n+1 


a 


(i,P) 


1 + 


+ vP + 


9e 


(*>p) 


7 


E 

p&P 


E, 


(q,p),{i,p) \ x n 


(<hP) — X (^P. 




■Sc E [ A «JU‘,r) H ( x n J) - - V *)] + A., 


(3) 


(d,i)€S 

(/,p)eP 


(i,p) 

Vn+1 
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where S (i — 1,2,... ,S ) is the total number of neurons in each small-world, 
P (p = 1, 2,..., P) is the number of cortical areas (P = 65 according to Fig. 
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[2D, g e is the electrical coupling strength, g c is the chemical coupling strength, 
cd*-p) i s the non-linearity parameter of the Rulkov map with values randomly 
coupled in the interval [4.1,4.4], a = 0.001, and p = —1.25. The term A n 
is an external excitatory perturbation which activates spikes in randomly 
chosen neurons. The adjacency matrices -E( g ,p),(i, p ) and A(d,f),(i,p) are the elec¬ 
tric and chemical connections, respectively. They have elements with value 
equal to 1 when neuron (■ q,p ) connects electrically with neuron (i,p), and 
neuron (rf,/) connects chemically with neuron (i,p). H(x) is the Heaviside 
step function, where 9 = —1.0 is the presynaptic threshold for the chemical 
synapse. When the presynaptic neuron voltage is above 9 , the postsynaptic 
neuron receives an input. This way, 9 is related to the sharp voltage response 
of the presynaptic terminals. The constant V s denotes the reversal potential 
associated with the synapse, that is defined by the nature of the postsynaptic 
ionic channels. The synapse will be excitatory if V s is higher and inhibitory 
if V s is lower than a specific range 35[. Regarding the network of Rulkov 


neurons, for excitatory synapses V s = 1.0 and for inhibitory V s = —2.0. 


3. Burst synchronisation without external perturbation 

Bursts of spikes and oscillatory patterns of neuronal activity have been 
observed in the central nervous system, ranging from slow to fast oscillations 


36]. In fact, synchronised bursting has been verified in EEG recording of 


electrical brain activity. A type of synchronisation is the burst phase syn¬ 
chronisation [371 ]. 

Burst phase synchronisation is studied through the definition of a phase 
for neuronal bursting, defined by the slow variable. We consider that a burst 
begins when the slow variable y n , shown in Fig. QJb), has a local maximum 
happening in a time n*,. The duration of the burst, Uk+i — Uk, depends on the 
variable x n and fluctuates in an irregular behaviour as long as x n undergoes 
irregular evolution. Then, we define a phase describing the time evolution 
within each burst, varying from 0 to 27T as n evolves from to rik+i, 


= 2tt k + 2it- 


n - n k 
^fc+1 W'k 


(5) 


where k is an integer. 

We use the Kuramoto’s order parameter as a diagnostic of the burst phase 
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synchronisation, that is defined as 

4> = H™ exp(it®) = A Y, exp(i0“’ , ' ) ), (6) 

where R n is the amplitude and <f> n is the angle of a centroid phase vector for 
an one-dimensional network with periodic boundary conditions. // represents 
one of the four cognitive areas, l = 1 for visual, / = 2 for auditory, l = 3 
for somatosensory-motor, and l = 4 for frontolimbic. IV) is the number of 
neurons of each area. represents the phase of the neurons j belonging 
to the c cortical area If the bursting phases (pn' 1 ^ are uncorrelated, the 
summation in Eq. (Jf)J) is small and Rn -C 1. Whereas R$ = 1 when the 
cortical network area is in a completely burst phase synchronised state. 

We are interested in studying the role of the electrical and chemical cou¬ 
pling strength at the level of bursting synchronisation. For this reason, we 
use the time averaged order parameter magnitude, given by 

(?) 

72—1 


where T is the time interval, as a measure of synchronisation in the network. 
If the bursting dynamics in the area 7/ is globally synchronised, we obtain 



Figure [3] shows the time averaged order parameter in colour scale as a 
function of the electrical and chemical coupling strength. Due to the fact 
that the neurons are not identical the neuronal network does not present 
a completely phase synchronised state (R^ ~ 1). However, the network 
exhibits strong synchronisation for R (l> > 0.9, corresponding to the white 
region. All the cognitive areas, except for the auditory, present strong syn¬ 
chronisation (white region). 

The auditory region does not present high levels of bursting synchro¬ 
nisation. This effect is due to the complex network topology, where the 
synchronisation patterns are controlled by input intensities among neurons, 
and also among the corticocortical areas |15j. To show that we define the 
mean held of the auditory area 
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Figure 3: (Colour online) Average order parameter of the cat’s cognitive brain regions. 
We consider 50000 iterations where 20000 were transient. 


and the mean field of the N OVLt inputs coming from other areas 



cdO 

On 


N t 


out 


where 


(9) 


3? =E [A< dMi:r) H(x^-6)(x^-V,)], (10) 

(i,p) er 2 
(dj)£h,h,l3 

where N out represents the number of non null elements in the matrix A^,f),(i,p) 
for (i,p) G / 2 and (d, f) G I\,l2,h- We assume that each neuron (i,p) 
belonging to a set / 2 in the auditory area is connected to a neuron (d. /) in 
a set L belonging to another cognitive area. Due to spatial distance between 
areas |38| we consider that the inputs from other areas are through chemical 
synapses [39 : ]. 

Aiming to demonstrate that the absence of high levels of synchrony in 
the auditory area is an effect of the network, we compare the mean fields 
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Figure 4: (Colour online) (a) Dynamics of the auditory fast variable mean field. Evolution 
of the external mean field on auditory area from (b) visual, (c) somatosensory-motor, and 
(d) frontolimbic areas. We consider g e = 0.05 and g c = 0.015. 

with Cn \ Cn \ and Cn\ Figure |4j a ) exhibits the time evolution of 
the mean field Af ?l ' for N 2 = 1000. The irregular time evolution, that it is 
associated with a non synchronised behaviour, is influenced by the stimuli 
from the visual (Fig. 0b), the somatosensory-motor (Fig. 0b), and the fron¬ 
tolimbic (Fig. 0Ji) areas. This figure shows that the collective behaviour of 
the auditory area is correlated with the input from the frontolimbic. If the 
auditory area is isolated from other areas, by making A( d j)eh,i 3 ,u-,{i,p) = 0 it 
will present global bursting synchronisation with a frequency around 0.0027. 
The frequency is obtained from bursts. Connecting the auditory area with 
the rest causes a suppression of bursting synchronisation, as well as the ap¬ 
pearance of one more frequency with value equal to 0.0025 (Fig. 0b). We can 
see through Figure 01(b) that the new frequency appears by a resonant effect 
caused by the oscillations of the stimulus from others areas. The visual and 
somato-motor areas present the frequency equal to 0.0025 with magnitude 
equal to 0.0005 and 0.001, respectively. The frontolimbic area also exhibits 
this value of frequency, but the magnitude is larger than visual and somato¬ 
motor areas. Therefore, the low levels of synchrony in the auditory area is 
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mainly due to the new frequency of the input from frontolimbic area. 



0.005 
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Figure 5: (Colour online) Frequencies of the (a) auditory fast variable mean field, and (b) 
external mean field on auditory area, considering g e = 0.05 and g c = 0.015. 


4. External perturbation inducing synchronisation 


It has been experimentally found that the incidence of light with a suit¬ 
ably chosen frequency on a group of neurons is able to alter their spiking 


activity. This phenomenon was verified by Boyden and collaborators [2a, 27 


by means of blue light they activated spikes in neurons genetically modified, 
and also suppressed spikes through yellow light. With this in mind, we con¬ 
sider an external perturbation (A n ) acting on 100 neurons randomly chosen 
in the same cognitive area. This perturbation, aiming to simulate blue light 
optical stimulus, has excitatory effect on the perturbed neurons making them 
spike independent of their previous state. 

We apply an external perturbation on the visual area to verify the syn¬ 
chronisation behaviour. Figure [6] shows the effects of this perturbation in 
the cognitive areas. Synchronisation is fully suppressed in the visual area. 
The somatosensory-motor and frontolimbic areas do not present significant 
alterations due to the perturbation in the visual area. This behaviour is due 
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Figure 6: (Colour online) Order parameter of the cat’s cognitive brain regions with an 
external perturbation on the visual area. We consider 50000 iterations, and 20000 transient 
iterations. 


to the small number of connections between these areas and the visual area. 
On the other hand, the auditory area exhibits a change in its synchronous 
behaviour. Without an external perturbation the auditory area does not 
present high levels of synchronisation, shown in Figure |3l but with the per¬ 
turbation in the visual area it is possible to observe strong synchronisation 
domains in the auditory area. 

In order to understand the perturbation effect in the auditory area we 
calculate the time evolution of the mean held , and the mean held from 
inputs coming from others areas Cn\ Cn \ and Cn^ into the area / 2 . Figure 
[Ha) shows that the mean held of the auditory area has a regular behaviour, 
whereas with no external perturbation the behaviour is irregular. Comparing 
the results of Figure |7] for the external inputs of the mean held on auditory 
from visual (b), somatosensory-motor (c), and frontolimbic (d) with the re¬ 
sults of Figure [4] we can see that the inputs from the visual area becomes 
frequency-locked with the mean held MA ' of the auditory area. MA and 
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Cn present both the same frequency, approximately 0.0028, as it can be 
seen in Fig. [8] 



n 

Figure 7: (Colour online) (a) Dynamics of the auditory fast variable mean field. Evolution 
of the external inputs mean field on the auditory area from (b) visual, (c) somatosensory- 
motor and (d) frontolimbic areas with an external perturbation on the visual area, con¬ 
sidering g e = 0.05 and g c = 0.015. 

Perturbations applied in auditory area do not alter significantly the syn¬ 
chronous behaviour of the other areas. Moreover, we have verified which 
perturbations in the somatosensory-motor and frontolimbic areas produce 
the same effect on the network such as when they are applied in the visual 
area. 

5. Conclusion 

In this paper we studied burst phase synchronisation in a neuronal net¬ 
work with a topology according to the corticocortical connections of the cat 
cerebral cortex. When no perturbation is applied we verified that only the 
auditory area does not present synchronisation in the parameter space. This 
happens because of the large influence of the other areas. The ratio between 
inter connections and intra connections is less for the auditory area than 
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Figure 8: (Colour online) Frequencies of the (a) dynamics of the auditory fast variable 
mean field and (b) inputs mean field on auditory area for the system with external per¬ 
turbation on visual area. 


other areas. As the consequence the frontoiimbic area induces a frequency 
in the auditory area which helps to suppress burst phase synchronisation. 

We verified the suppression of burst phase synchronisation when a stimu¬ 
lus is applied on a cognitive area. However, the suppression in a specific area 
can affect another area. Our results have showed that an external pertur¬ 
bation applied into the visual area suppresses synchronisation in the visual 
area, but surprisingly it induces synchronisation in the auditory area. The 
same happens for perturbation applied in the somatosensory-motor and fron- 
tolimbic which induces synchronisation in the auditory area. 
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